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ABSTRACT 

We use cosmological hydrodynamical zoom-in simulations with the SPH code GASO¬ 
LINE of four haloes of mass M 200 ~ 10 13 Mq to study the response of the dark matter 
to elliptical galaxy formation. Our simulations include metallicity dependent gas cool¬ 
ing, star formation, and feedback from massive stars and supernovae, but not active 
galactic nuclei (AGN). At z = 2 the progenitor galaxies have stellar to halo mass ratios 
consistent with halo abundance matching, assuming a Salpeter initial mass function. 
However by z = 0 the standard runs suffer from the well known overcooling problem, 
overpredicting the stellar masses by a factor of > 4. To mimic a suppressive halo 
quenching scenario, in our forced quenching (FQ) simulations, cooling and star forma¬ 
tion are switched off at z = 2. The resulting z = 0 galaxies have stellar masses, sizes 
and circular velocities close to what is observed. Relative to the control simulations, 
the dark matter haloes in the FQ simulations have contracted, with central dark mat¬ 
ter density slopes d\ogp/d\og r ~ —1.5, showing that dry merging alone is unable to 
fully reverse the contraction that occurs at z > 2. Simulations in the literature with 
AGN feedback however, have found expansion or no net change in the dark matter 
halo. Thus the response of the dark matter halo to galaxy formation may provide a 
new test to distinguish between ejective and suppressive quenching mechanisms. 

Key words: cosmology: theory - dark matter - galaxies: elliptical and lenticular, cD 
- galaxies: formation - galaxies: evolution - methods: numerical 


1 INTRODUCTION 

Dissipationless simulations run in the concordance ACDM 
cosmology make robust predictions for the structure of cold 
dark matter haloes (e.g., Navarro et al. 1997, 2010; Bul¬ 
lock et al. 2001; Diemand et al. 2007; Maccio et al. 2007; 
Stadel et al. 2009; Zhao et al. 2009; Klypin et al. 2011; 
Dutton & Maccio 2014). One of these predictions is that 
CDM haloes should have “cuspy” central density profiles 
that scale as p(r) oc r~ a with a tn —1.2. Observations 
have yet to unambiguously find these cusps. Furthermore, 
a variety of studies based on gas dynamics, stellar dynamics 
and gravitational lensing on scales from dwarf galaxies to 
galaxy clusters often favor cored (a = 0) or shallow cusps 
(i.e., a < 0.5) instead (e.g., de Blok et al. 2001; Swaters 
et al. 2003; Sand et al. 2004; Goerdt et al. 2006; Kuzio de 
Naray et al. 2008; Walker & Penarrubia 2011; Oh et al. 2011; 
Newman et al. 2011; 2013). 

This lack of observational evidence for cuspy dark mat- 
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ter density profiles is often used as evidence against the 
CDM paradigm. However, a major complication in the com¬ 
parison between observations of dark matter density profiles 
and those predicted by cosmological N-body simulations is 
the uncertain impact of baryons. Baryonic processes can 
both increase and decrease the density profiles of dark mat¬ 
ter haloes. This currently prevents the observed structure of 
dark matter haloes from being used as a robust test of the 
cold dark matter model. A fully predictive theory for the 
structure of CDM haloes must take into account the effects 
of galaxy formation. 

If the accretion of baryons onto the central galaxy is 
smooth and slow then dark matter haloes should contract 
adiabatically (Blumenthal et al. 1986). This process can in¬ 
crease the density of dark matter haloes by an order of mag¬ 
nitude. Other processes can cause the dark matter halo to 
expand, such as transfer of energy/angular momentum from 
baryons to the dark matter via dynamical friction due to 
minor mergers (El-Zant et al. 2001, 2004; Nipoti et al. 2004; 
Jardel & Sellwood 2009; Johansson et al. 2009; Lackner & 
Ostriker 2010; Cole et al. 2011; Laporte & White 2015), or 
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galactic bars (Weinberg & Katz 2002; Sellwood (2008), but 
see McMillan & Dehnen 2005), and rapid mass loss/time 
variability of the potential due to supernovae (SN) / stellar 
/ active galactic nuclei (AGN) feedback (Navarro et al. 1996; 
Read & Gilmore 2005; Mashchenko et al. 2006, 2008; Peirani 
et al. 2008; Governato et al. 2010; Pontzen & Governato 
2012; Maccio et al. 2012; Martizzi et al. 2012, 2013). Each 
of these processes is likely to occur during galaxy forma¬ 
tion, but at present it is not clear which process (if any) 
dominates, and on which galaxy mass scales. 

Turning the problem around, if one can observationally 
measure the halo response, then in the context of ACDM, 
this gives clues to the dominant galaxy formation mecha¬ 
nisms. In particular, as we discuss in this paper, the halo 
response may provide a test for different galaxy quenching 
models. While it is known that most of the stars in mas¬ 
sive galaxies formed at high redshifts a > 2 (e.g., Thomas 
et al. 2005), the mechanism(s) responsible for shutting down 
star formation (quenching), and keeping it off (quiescence) 
are a subject of much debate. A popular idea is that feed¬ 
back of energy from the growth of supermassive black holes 
(SMBH) ejects cold gas from galaxies and heats the halo 
gas, preventing it from cooling onto the central galaxy. We 
refer to this class of models as “AGN quenching”. It is clear 
there is enough energy released from the AGN to quench 
star formation, of issue is whether the energy released by 
the AGN can couple efficiently to the surrounding gas (e.g., 
Cielo et al. 2014). 

An alternative idea is that cooling becomes very ineffi¬ 
cient, and effectively shuts down, once the halo mass reaches 
a critical mass of ~ 10 12 Mq. We refer to this class of mod¬ 
els as “halo quenching”. In practice the boundary between 
these two quenching mechanisms is blurry because AGN can 
also act as a heat source for the halo gas, although there are 
non-AGN heating sources such as young and old stellar pop¬ 
ulations (Kannan et al. 2014, 2015; Conroy et al. 2015). The 
key difference in these mechanisms is that AGN quenching 
involves the episodic (and possibly violent) removal of gas 
from the galaxy center, whereas halo quenching involves the 
suppression of gas cooling from the halo. 

Galaxy formation models implementing both mecha¬ 
nisms can broadly reproduce the colors and mass functions 
of present day galaxies (e.g., Croton et al. 2006; Cattaneo 
et al. 2006; Schaye et al. 2015). However, there is at present 
no clear way to observationally distinguish between these 
models. The correlation between SMBH mass and galaxy 
properties is often cited as evidence for a physical link be¬ 
tween SMBHs and host galaxies, although (Jahnke & Maccio 
2011) showed that such a correlation can naturally arise out 
of hierarchical merging. 

As discussed above, in terms of formation mechanisms, 
there is one process that makes haloes contract: dissipative 
gas accretion (Blumenthal et al. 1986), which must occur 
to form an elliptical galaxy. Two processes have been inves¬ 
tigated that make haloes expand: mass outflows driven by 
AGN feedback (Martizzi et al. 2012); and dynamical friction 
between baryons and dark matter during the galaxy assem¬ 
bly process (El-Zant et al. 2001). Dry mergers are likely to 
be more important at late times, while AGN feedback is 
likely to be more important at early times. 

The goal of this study is to determine the impact of 
dissipationless galaxy assembly on the structure of ACDM 


Table 1. Cosmological parameters. Column (1) cosmology ID. 
Column (2), fl m , present day matter density. Column (3), C,\, 
dark energy density. Column (4), , baryon density. Column (5), 

Ho, Hubble parameter. Column (6), erg, power spectrum normal¬ 
ization. Column (7), n, power spectrum slope. 


Cosmology f2 m 1 Ho erg n 

(1) (2) (3) (4) (5) (6) (7) 

WMAP5 0.2580 0.7420 0.0438 0.720 0.796 0.963 

WMAP7 0.2748 0.7252 0.0458 0.702 0.816 0.968 


haloes. In particular, we wish to know whether dry merg¬ 
ing can reverse the effects of halo contraction that are ex¬ 
pected to occur in the high redshift progenitors. To achieve 
this we use fully cosmological hydrodynamical simulations of 
the formation of M 200 ~ 10 13 Mq haloes using the smoothed 
particle hydrodynamics (SPH) code GASOLINE. We use the 
MaGICC (Stinson et al. 2013) star formation and stellar 
feedback model which has been shown to produce realistic 
disk galaxies (Brook et al. 2012) and successfully matches 
the observed galaxy formation efficiencies across cosmic time 
in haloes less massive than ~ 10 12 Mq (Stinson et al. 2013; 
Wang et al. 2015). Our study improves on previous analytic 
and N-body calculations (e.g., El-Zant et al. 2001; Lackner 
& Ostriker 2010) in that our mass assembly histories are 
fully cosmological, and compared to previous cosmological 
hydro simulations (e.g., Johansson et al. 2009) our progen¬ 
itor galaxies have realistic sizes and galaxy formation effi¬ 
ciencies. 

This paper is organized as follows: the simulations in¬ 
cluding sample selection, hydrodynamics, star formation 
and feedback models are described in fj2j Results relating to 
global parameters such as stellar masses, halo masses, galaxy 
sizes and circular velocities are presented in fj3j The radial 
mass profiles including the response of the dark matter to 
galaxy formation are presented in ||4] We discuss implica¬ 
tions for our results in an and give a summary in ij6] 


2 SIMULATIONS 

The simulations presented here are fully cosmological 
“zoom-in” simulations of galaxy formation run in a flat 
ACDM cosmology. Cosmological parameters (see Table [T]l 
are based on the Wilkinson Microwave Anisotropy Probe 
(WMAP) 5th year (Komatsu et al. 2009) and 7th year (Ko¬ 
matsu et al. 2011) results. Haloes are selected from two 
parent dark matter only simulations run with the pkd- 
GRAV tree-code (Stadel 2001): a 90 Mpc box from Maccio 
et al. (2008) using the WMAP5 cosmology, and a 80 Mpc /h 
box from Penzo et al. (2014) using the WMAP7 cosmology. 

2.1 Sample selection and initial conditions 

We select four haloes to re-simulate at higher resolutions. 
The two selection criteria are (1) a present day halo mass of 
~ 10 13 Mq; and (2) no large structures present within three 
virial radii. Table [2]lists the dark matter and baryon particle 
masses and force softenings (in comoving kpc). To test for 
numerical convergence, one of the haloes (halo4) has been 
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Table 2. Simulation parameters. Column (1), name of initial conditions. Column (2), M 200 , present day halo mass from adiabatic gas 
simulation. Column (3), dark matter particle mass. Column (4), m gas , initial gas particle mass. Column (5), edarki dark matter 

particle force softening in comoving kpc. Column (6), e ga s, gas (and star) particle force softening in comoving kpc. Column (7), iVdarki 
number of high-resolution dark matter particles. Column (8), lV gas , initial number of gas particles. Column (9), n t h, threshold for star 
formation. Column (10) name of cosmology - see Tablejl] Column (11) symbol used in figures. 


Name 

(i) 

logic M 200 
[M 0 ] 

(2) 

logio Uplajk 

[M 0 ] 

(3) 

logio ™gas 

[M 0 ] 

(4) 

^dark 

[ckpc] 

(5) 

e gas 

[ckpc] 

(6) 

Vla.rk 

Million 

(7) 

Ngas 

Million 

(8) 

™th 

[cm -3 ] 

(9) 

Cosmology 

(10) 

symbol 

(ii) 

halol 

13.15 

6.91 

6.23 

1.61 

1.61 

2.49 

2.49 

9.6 

WMAP5 

triangle 

halo2 

13.19 

6.91 

6.23 

1.61 

1.61 

2.49 

2.49 

9.6 

WMAP5 

square 

halo3 

13.20 

6.91 

6.23 

1.61 

1.61 

2.49 

2.49 

9.6 

WMAP5 

pentagon 

halo4.0 

13.36 

8.13 

7.43 

4.06 

1.81 

0.34 

0.34 

1.16 

WMAP7 

circle 

halo4.1 

13.35 

7.23 

6.53 

2.04 

0.910 

2.68 

2.68 

1.16 

WMAP7 

circle 

halo4.2 

13.33 

6.70 

6.00 

1.36 

0.606 

9.06 

9.06 

1.16 

WMAP7 

circle 

halo4.3 

13.42 

6.32 

5.63 

1.02 

0.455 

21.5 

21.5 

1.16 

WMAP7 

circle 


run at four different resolution levels, with particle masses 
varying by a factor of 64, and force softenings by a factor 
of 4. Note that for halo4 the slightly higher halo mass in 
the highest resolution run (halo4.3) is due to a merger event 
that occurs slightly earlier in this simulation compared to 
the lower resolution runs. 

Fig. Q] shows the dark matter mass resolution of our 
simulations (red filled symbols) compared to a number of 
state-of-the-art simulations in the literature of comparable 
mass haloes. We only consider simulations that have been 
run all the way to redshift z — 0. We note that the simula¬ 
tions from Feldmann & Mayer (2015) are higher resolution, 
with a dark matter particle mass of 7.9 x 10 5 Mq, however, 
these simulations are only run to z = 2, where the halo 
mass is ~ 3 x 1O 12 M 0 and thus they do not appear in our 
comparison. 

The highest resolution simulation of a ~ 1O 13 M 0 halo 
in the literature is from the 25 Mpc box from the EAGLES 
project (Schaye et al. 2015, open circles). This simulation 
has a dark matter particle mass of ttidm = 1.21 x 1O 6 M 0 
and thus ~ 7 million particles for a ~ 1O 13 M 0 halo. Compa¬ 
rable resolution can be obtained with zoom-in simulations 
with a fraction of the computational cost. The FIRE project 
(Hopkins et al. 2014, cyan square) includes a simulation with 
a dark matter particle mass of moM = 2.26 x 1O 6 M0, and 
roughly 5 million dark matter particles within the virial ra¬ 
dius. For comparison our highest resolution runs of halo4 
have 3.6 and 10.6 million dark matter particles inside R 200 , 
and thus are among the highest resolution elliptical galaxy 
simulations performed to date. 

At the fiducial resolution our simulations have ~ 1.5 
million dark matter particles inside Raoo- This is similar to 
that of the high resolution halo in Feldmann et al. (2010, 
blue pentagon), the haloes in Dubois et al. (2013, magenta 
hexagons), together with the large volume (~ 100 Mpc) 
simulations by the ILLUSTRIS (Vogelsberger et al. 2014) 
and EAGLE collaborations (dashed lines). Earlier “zoom- 
in” simulations from Oser et al. (2010) and Feldmann 
et al. (2010) have particle masses of ~ 3.6 x 1O 6 M 0 , and 
thus ~ 300, 000 dark matter particles per halo. 
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Figure 1. Dark matter mass resolution of our simulations (red 
filled symbols) compared to state-of-the-art simulations of ellip¬ 
tical galaxies in the literature (open symbols). 


2.2 Hydrodynamics 

Our standard simulations use the same baryonic physics that 
was used in the MaGICC project (see Stinson et al. 2013), 
based on the smoothed particle hydrodynamics (SPH) code 
GASOLINE (Wadsley et al. 2004). 

Cooling via hydrogen, helium, and various metal-lines 
in a uniform ultraviolet ionizing background is included as 
described in Shen et al. (2010) and was calculated using 
cloudy (version 07.02; Ferland et al. 1998). These calcula¬ 
tions include photo ionization and heating from the Haardt 
& Madau (2005) UV background and Compton cooling in 
the temperature range 10 to 10 9 K. In the dense, interstel¬ 
lar medium gas, we do not impose any shielding from the 
extragalactic UV field as the extragalactic field is a reason- 
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Figure 2. SUNRISE images of simulations at redshift z = 0. Each box is 50 kpc on a side. Upper two panels show results for standard 
(overcooled) simulations, while lower two panels show forced quenching simulations. From left to right: halol, halo2, halo3, and halo4.2. 


able approximation in the interstellar medium. Diffusion of 
metals and thermal energy between particles has been im¬ 
plemented as described in Wadsley et al. (2008). 

Stars form from cool dense gas that has reached a tem¬ 
perature of T < 1.5 x 10 4 K and a density of n > n t h- 
For halol-3 we adopt n t h = 9.6 cm -1 following MaG- 
ICC (Stinson et al. 2013). For halo4 we adopt nth = 1-16 
cm -1 , which is a conservative estimate of the maximum 
gas density that can be resolved and is calculated using 
n t h = 32(m gaa /5)/e| aa , where m gas is the initial gas par¬ 
ticle mass, m gaa /5 is the minimum gas particle mass, and 
e gaa is the gas force softening. The factor 32 is the num¬ 
ber of SPH smoothing elements. This difference in star for¬ 
mation threshold has no noticeable impact on the struc¬ 
ture of the galaxies or dark matter halo. Star formation fol¬ 
lows the Kennicutt-Schmidt law with 10 percent efficiency 
of turning gas into stars during one dynamical time (Stinson 
et al. 2006). The stellar mass distribution in each star parti¬ 
cle follows the Chabrier (2003) stellar initial mass function 
(IMF). 

Massive stars explode as Type II SN and deposit an 
energy of Esn = 10 51 erg into the surrounding gas. Cooling 
for gas particles subject to SN feedback is delayed based on 
the subgrid approximation of a blast wave as described in 
Stinson et al. (2006). Furthermore, radiation energy from 
massive stars is considered since molecular clouds are dis¬ 
rupted before the first SN explosion (which happens after 4 
Myr from the formation of the stellar population). We as¬ 
sume that 10 percent of the total radiation energy is coupled 
with the surrounding gas. The inclusion of this early stellar 
feedback (ESF) reduces star formation before SN start ex¬ 
ploding. Thus, after the ESF heats the gas to T > 10 6 K, 
the gas rapidly cools to 10 4 K, which creates a lower density 
medium than if the gas were allowed to continue cooling un¬ 
til SN exploded. Stinson et al. (2013) and Wang et al. (2015) 


show how this feedback mechanism limits star formation to 
the amount prescribed by halo abundance matching (e.g., 
Behroozi et al. 2013) at all redshifts for haloes of present 
day virial mass ~ 10 10 — 10 12 Mq. 

For each initial condition we run three simulations that 
differ only in their treatment of cooling and star formation. 

• Standard: cooling, star formation, and stellar feedback 
at all times following MaGICC (Stinson et al. 2013). Feed¬ 
back from AGN is not included. 

• Forced Quenching (FQ): same as standard down to 
z ~ 2.1, after which the simulation is evolved with adiabatic 
gas until z = 0 - i.e., there is no further star formation. 

• No Cooling: Gas is adiabatic at all times, and thus 
does not become cool and dense enough to form stars. These 
simulations act as a control for the effects of galaxy forma¬ 
tion (cooling, star formation, and feedback) on the dark halo 
structure. 

The forced quenching simulations act as a limiting case 
where cooling is shut down and star formation stops rapidly. 
It allows us to study the effects of dry merging on the struc¬ 
ture of the galaxy and dark matter halo with cosmologi- 
cally consistent initial conditions and merger histories. As 
we show below this simple prescription allows one to re¬ 
produce the stellar mass vs halo mass relation from abun¬ 
dance matching, which is not possible when adopting stan¬ 
dard metal line cooling, and no additional heating sources 
(such as AGN feedback). We note that previous simula¬ 
tions that come close to matching the stellar mass vs halo 
mass relation without AGN feedback (Oser et al. 2010; Feld- 
mann et al. 2010; Feldmann & Mayer 2015) do so because 
they adopt primordial metallicity gas in the calculation of 
the cooling rate. Including realistic halo gas metallicities in¬ 
creases the cooling rate by a factor of ~ 2 — 3 (e.g., Dutton 
& van den Bosch 2012), and thus would result in overcooled 
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galaxies. It has been shown that ionizing radiation from 
young and old stellar populations can significantly reduce 
cooling rates in 1O 12 M0 haloes (Kannan et al. 2014). Addi¬ 
tionally Conroy et al. (2015) argue that heating provided by 
the winds of dying low-mass stars is capable of suppressing £7 

cooling of hot gas for a Hubble time in halo masses above 2 

~ 10 12,5 Mq at redshifts below a ~ 2. Thus, there is moti- 
vation for studying the impact of inefficient gas cooling on 
elliptical galaxy formation. ^ 

cn 

O 

2.3 Derived galaxy and halo parameters in 

Haloes in our zoom-in simulations were identified using 
the MPI+OpenMP hybrid halo finder AHI0 (Knollmann & 

Knebe 2009; Gill et al. 2004). AHF locates local over-densities 
in an adaptively smoothed density field as prospective halo 
centers. The virial masses of the haloes are defined as the 
masses within a sphere containing A = 200 times the cosmic 
critical matter density. The virial mass, radius, and circular 
velocity are denoted M 200 , R200, and V 2 oo- At 2 = 0 there 
is, by construction, one central halo in the zoom-in region. 

However, at 2 ~ 2 there are several well resolved progeni¬ 
tor galaxies which we also consider in the global parameter 
evolution plots. 

The mass in stars, M s t ar , is measured within a sphere 
of radius, r ga 1 = 0.2R2OO- The stellar half-mass radius, 77/2, 
encloses within a sphere half of the stellar mass within 
r ga i. The half-mass circular velocity is defined as V\/ 2 = 
yjGM(r 1 / 2 )/ri/ 2 , where M(ri/ 2 ) is the total mass within a 
sphere of radius r 1 / 2 . The (mass weighted) slope of the total 
mass density profile, 7', and the “inner” slope of the dark 
matter density profile, a, are measured between 0.01 and 
0.02 R 2 oo- The choice of this scale is motivated by the follow¬ 
ing: it is used by previous studies (e.g., Di Cintio et al. 2014); 
it is resolved in our simulations; and it corresponds to the 
average half-light radii of galaxies (Kravtsov 2013) and is 
thus observationally accessible. 

Images of our simulated galaxies at a = 0 are shown 
in Figs. [2] Each image is 50 kpc on a side and was created 
using the Monte Carlo radiative transfer code SUNRISE (Jon- 
sson 2006). The image brightness and contrast are scaled 
using arcsinh as described in Lupton et al. (2004). All simu¬ 
lated galaxies have red colors and elliptical like morphology. 

Some show signs of recent mergers in the form of shells and 
streams. The main apparent difference between the standard 
(upper panels) and forced quenching (lower panels) simula¬ 
tions is a reduction in surface brightness (the same scale is 
used in all images). As is visible, and we show quantitatively 
below, the half-light sizes are roughly the same in each type 
of simulation (standard vs FQ). 


3 GLOBAL PROPERTIES 

Before we discuss the response of dark matter haloes to 
galaxy formation in § [I] we present several global proper¬ 
ties of our simulated galaxies and compare them to obser¬ 
vations. At 2 = 2.1, the standard simulations compare well 
with observed galaxies. However, once they are evolved to 

1 http://popia.ft.uam.es/AMIGA 
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Figure 3. Stellar mass vs halo mass relation for simulations at 
redshift 2 = 0 and the progenitor galaxies at 2 = 2.1. Open sym¬ 
bols show standard simulations (with cooling and star formation 
at all redshifts) at 2 = 0 (black) and 2 = 2.1 (blue for most mas¬ 
sive progenitors, cyan for other progenitors). Solid red symbols 
show forced quenching simulations (which have no cooling or star 
formation since 2 = 2.1). The dashed lines show the abundance 
matching relations from Behroozi et al. (2013) at redshifts 2 = 0 
(red) and 2 = 2.1 (blue), the solid lines show the corresponding 
relations for stellar masses higher by 0.23 dex, corresponding to 
a Salpeter IMF. The shaded ares shows the 1 a scatter in stellar 
mass at fixed halo mass. 

2 = 0, the standard simulations have too many stars, are 
too dense, and have mass profiles that are too centrally con¬ 
centrated. Thus, they are poor candidates for studying halo 
contraction. In contrast, we show that the forced quenching 
simulations share many properties with observed galaxies at 
2 = 0, so it is more interesting to study them as templates 
for elliptical galaxy formation. 


3.1 Stellar mass vs halo mass 

One of the most fundamental properties of a galaxy or dark 
matter halo is its mass. The relation between the mass in 
stars, Mstar, and the virial mass, M200, of our simulated 
galaxies is shown in Fig. [3] Results are shown at redshift 
2 = 0 for standard (black open symbols) and forced quench¬ 
ing (red filled symbols) and at 2 = 2.1 (blue open symbols 
show the most massive progenitor, while cyan open symbols 
show other progenitors containing more than 20,000 total 
particles). The lines show the relations from halo abundance 
matching from Behroozi et al. (2013), with the shaded re¬ 
gion showing the ler intrinsic scatter. Two assumptions for 
the observed stellar masses are shown: A Chabrier (2003) 
initial mass function (IMF, dashed) as found in the Milky 
Way; and Chabrier+0.23 dex (solid), which corresponds to 
a Salpeter (1955) IMF and is likely a better approximation 
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Figure 4. Scaling relations between galaxy size, rq/ 2 , circular velocity at the half-mass size V^/ 2 , and galaxy stellar masses, M s t ar , 
at z = 2.1 and 2 = 0. The observed relations for quiescent galaxies are shown at z = 0 (red/pink) and z = 2.1 (blue/grey) for both 
Chabrier (dashed) and Salpeter (solid) IMFs. The shaded regions show the la scatter. At 2 = 0 the standard simulations are shown with 
black open symbols, and the forced quenching simulations are shown with red filled symbols. At z = 2.1 the most massive progenitors 
are shown with blue symbols, and the lower mass progenitors with cyan. The cyan lines show the observed size-mass and velocity-mass 
relations for star forming galaxies at z = 2.1. 


for the centers of massive elliptical galaxies (e.g., Conroy & 
van Dokkum 2012; Dutton et al. 2013b). 

At redshift 2 = 0 the standard simulations over predict 
the Chabrier IMF stellar masses by an order of magnitude 
— this is the classic overcooling problem. Even if one adopts 
a Salpeter IMF, the stellar masses are still over predicted by 
a factor of ~ 4. Given the 1 <t intrinsic scatter in the stellar 
masses at fixed halo masses is ~ 0.2 dex (More et al. 2011; 
Reddick et al. 2013), real galaxies offset from the median 
relation by a factor of ~ 4 will be extremely rare, if they 
even exist. Thus there is no escape from the conclusion that 
these simulations have formed too many stars. At 2 = 2.1, 
however, the standard simulations have only a factor of ~ 2 
too many stars compared to a Chabrier IMF. Thus when 
adopting an observed Salpeter IMF the simulations provide 
a good match to the halo abundance matching results. Note 
that while there is a formal inconsistency here between the 
IMFs used in our simulations (Chabrier) and observation 
(Salpeter), in practice the feedback efficiency parameters in 
the simulations can be simply rescaled. 

When we disable cooling and force quenching at red- 
shift 2 = 2.1 the resulting stellar masses at 2 = 0 are 
consistent with halo abundance matching constraints, pro¬ 
vided we also adopt a Salpeter-like IMF at 2 = 0. Thus, 
our simulations show that no new stars are needed to be 
formed since redshift 2 ~ 2 in haloes of present day mass of 
M 200 ~ 1O 13 M 0 . This result supports a two phase galaxy 
formation picture (e.g., Oser et al. 2010). At halo masses 
below M 200 ~ 10 12 Mq central galaxies form the majority of 
stars in situ, with the efficiency being regulated by stellar 
feedback. Above M 200 ~ 10 12 Mq in situ star formation is 


suppressed, and the growth in stellar mass is dominated by 
dissipationless mergers. 

A consequence of the dissipationless assembly is that 
whatever IMF characterizes the 2 = 0 galaxies must also 
describe their most massive progenitors. In our simulations 
progenitor galaxies with M s t ar > 10 10 Mq at 2 = 2.1 con¬ 
tribute 90% of the present day mass. Since the lower mass 
progenitors (M s t ar < 10 10 Mq) don’t contribute significantly 
to the 2 = 0 central galaxy mass, this self-consistency ar¬ 
gument does not constrain their IMF. The stars from the 
lower mass progenitors preferentially end up at large radii, 
which would result in a radial IMF gradient if the low mass 
galaxies have normal IMFs. 

3.2 Galaxy sizes and circular velocities 

The next global properties we consider are the half-mass size 
of the stellar mass distribution, 7q/2, and the circular veloc¬ 
ity at ri/ 2 . Fig. [4] shows the size-mass (left) and velocity- 
mass (right) relations at 2 ~ 0 and 2 ~ 2. 

For the size-mass relation we compare to observations 
from van der Wei et al. (2014) — who measure 2D major 
axis half-light radii in rest frame V-band light for both qui¬ 
escent and star forming galaxies. The observed relations for 
non-star-forming galaxies at 2 = 2.1 and 2 = 0 are shown 
with blue and red lines, respectively, with the shaded re¬ 
gion showing the la intrinsic scatter. As in Fig. [3] the solid 
lines assume a Salpeter IMF, while the dashed lines assume 
a Chabrier IMF. At 2 = 2.1 we also show the size-mass 
relation for star forming galaxies with cyan lines. 

In our simulations we measure spherical 3D half-stellar 
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Figure 5. Change in galaxy sizes, circular velocities, and stellar masses, between z = 2.1 and z = 0. of the most massive progenitor for 
the forced quenching (red, filled symbols) and standard (black, open symbols) simulations. The left panel shows the change in size vs 
change in mass, while the right panel shows the change in velocity vs change in mass. The red lines show the evolution required to stay 
on the size-mass and velocity-mass relations, with the shaded region showing the scatter in the 2 = 0 relation. For reference, the dotted 
lines show the approximate evolutionary paths for dry major and minor mergers. 


mass sizes. To correct for projection effects we multiply the 
observed sizes of ETGs by 4/3. No correction is applied to 
LTGs as the high disk fractions suggest a minimal difference 
between projected and spherical aperture sizes. We leave an 
investigation of projection effects, mass-to-light ratio varia¬ 
tions and axis ratios to a future study. The purpose of our 
current comparison is to determine if the simulated galaxies 
have reasonable structural properties. The observed relation 
between 3D half-light size, stellar mass (assuming a Chabrier 
IMF) and redshift is thus 


Rl/2 

[kpc] 


5.73 


M S i 


[5 x 10 10 Mq] 


H(z) 


(1) 


For the velocity-mass relation we use the 2 = 0 Faber & 
Jackson (1976) relation for quiescent galaxies from Gallazzi 
et al. (2006). We convert velocity dispersions into circular 
velocity using V) ; 2 = 1.5<r e (Dutton et al. 2011b; Cappel- 
lari et al. 2013). For the evolution we adopt the scaling of 
c(Mstar) oc (1 + 2 ) 0 ' 44 from the observed compilation of Oser 
et al. (2012). The observed relation between circular veloc¬ 
ity, stellar mass (assuming a Chabrier IMF) and redshift is 
thus 


Vl/2 

[kms -1 ] 


220 


M* 


[5 x 1O 1O M 0 ] 


0.29 

(1 + zf M 


( 2 ) 


For star forming galaxies galaxies at z ~ 2 we use the 
velocity-mass relation from Dutton et al. (2011a) using data 
from Cresci et al. (2009). 

The most massive galaxies at z = 2.1 in the standard 
simulations (blue symbols) have half-mass sizes of ~ 1 — 2 
kpc and circular velocities of ~ 400 — 500kms -1 , consis¬ 


tent with the size-mass and velocity-mass relations of qui¬ 
escent galaxies (dark grey shaded regions). The lower mass 
progenitors (cyan symbols) lie close to the corresponding 
relations for star forming galaxies (shown with cyan lines). 
Thus the standard simulations form a realistic population 
of progenitor galaxies, adding to the successes of MaGICC 
feedback model (Stinson et al. 2013) in reproducing global 
scaling relations of spiral galaxies (Brook et al. 2012; Wang 
et al. 2015). 

By 2 = 0 in both the standard and FQ simulations the 
most massive galaxies have grown substantially in size and 
mass. The standard simulations (black open symbols) have 
grown roughly along the 2 = 2.1 size-mass relation and are 
thus too dense compared to 2 = 0 observations. The FQ 
simulations fall parallel to the observed size-mass relation, 
being offset on average by factor of ~ 1.5. Observational ef¬ 
fects that may contribute to this discrepancy are discussed 
below. Here we note that numerical resolution may play a 
role, as the sizes of halo4 have not shown convergence when 
we increase the mass resolution by a factor of 8 (from levell 
to level3). By contrast, stellar masses are stable for simula¬ 
tions with more than a million particles. 

The changes since 2 = 2.1 in the sizes, circular veloc¬ 
ities, and stellar masses are more clearly shown in Fig. [5] 
For both standard and FQ simulations the sizes increase on 
average by ~ 0.6 dex, albeit with more variation in the FQ 
simulations. This is interesting given that the mass evolu¬ 
tion differs by a factor of ~ 4 between standard and FQ 
simulations. However, the velocity evolution is different: a 
decrease of ~ —0.1 dex for the FQ simulations, and an in¬ 
crease of ~ 0.15 dex for the standard simulations. 

The dotted lines in Fig. [5] show the changes expected 
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Figure 6. Logarithmic slope of the total mass profile between 
1 and 2 % of the virial radius, 7 ' vs the stellar surface density 
within the effective radius, E s t ar . Points with error bars show 
observations from the SLACS survey (Auger et al. 2010b). Stan¬ 
dard simulations are shown with black open symbols, while forced 
quenching simulations are shown with red solid symbols. 


for dry major (R oc M, V = const) and minor mergers 
(R oc M 2 , V oc A/ -1 / 2 ), (e.g., Bezanson et al. 2009; Hopkins 
et al. 2010). All but one simulation (halo2, squares) falls 
within these expected scalings. In halo2 the size increases 
by a factor of 2.5 with very little increase in stellar mass 
which points towards an adiabatic expansion scenario. 

The solid lines in Fig. [5] show the change in size and 
velocity required to remain on the size-mass and velocity- 
mass scaling relations between 2 = 2.1 and 2 = 0. For no 
mass change, sizes need to increase by 0.61 dex and veloci¬ 
ties need to decrease by —0.22 dex. If the stellar masses in¬ 
crease, even stronger evolution is required to remain on the 
size-mass relation. Our FQ simulations have the observed 
amount of evolution in velocities, but don’t have as much 
evolution in sizes as observed, they are still too small by a 
factor of ~ 1.5. There are several considerations to keep in 
mind when interpreting the small size: 

1) Population growth. If there is substantial growth in 
the number of quenched galaxies at a given mass, then the 
scaling relation will be dominated by the new arrivals. Thus 
earlier quenched galaxies could be significantly offset from 
the relation without violating the scatter. 

2) Progenitor size bias. Star forming galaxies are ob¬ 
served to be larger than quenched galaxies. If the quenching 
process does not change the galaxy size, then galaxies that 
are added to the quenched population will have, on aver¬ 
age, larger sizes than existing quenched galaxies of the same 
mass. Thus the earlier quenched galaxies will not need to 
grow as much in size as the population appears to. 

3) M/L variations. In the simulations we are measuring 
sizes in stellar mass, while observations are for optical stel¬ 



Figure 7. Relation between circular velocity at the half-mass 
radius, V 1 / 2 , and virial radius, V 200 at z = 0. Forced quench¬ 
ing simulations are shown with red solid points, and standard 
simulations with black open points. The shaded region shows ob¬ 
servations from Dutton et al. (2010). 

lar light. Galaxies are known to have smaller sizes in longer 
wavelengths (e.g., van der Wei et al. 2014). For quiescent 
galaxies the implied conversion from rest frame V-band to 
K-band (which presumably traces stellar mass) is on aver¬ 
age 0.15 dex. Star forming galaxies have stronger gradients 
at later times, which results in stronger evolution in half- 
light than half-mass sizes by 0.14 dex. A similar effect for 
quiescent galaxies would thus explain the entire discrepancy 
between our simulations and observations. In addition to 
color gradients, IMF gradients could plausibly make sizes 
appear larger in light than they are in mass. For example, 
if the most massive 2 = 2.1 progenitor formed with a uni¬ 
formly heavy IMF, and then grew by adding galaxies with 
lighter IMF to the outer parts. We note there are recent 
observational hints of such radial IMF variations in massive 
quiescent galaxies (Martfn-Navarro et al. 2015). 

3.3 Mass density slopes 

A further observational test of the realism of the simulated 
galaxies comes from the slope of the total mass profile, 7 ', 
which is measured near the galaxy half-light radius. Fig. [ 6 ] 
shows 7 ' vs stellar surface density, E star = M st a. r /(27rJ? 2 ) 
for our simulations and observations at 2 ~ 0. Points with 
error bars show observations using strong lensing from the 
SLACS survey (Auger et al. 2010b), where the stellar density 
is calculated assuming a Salpeter IMF. In our simulations 
we measure 7 ' between 1 % and 2% of the virial radius, as 
this corresponds to the typical half-light sizes of quiescent 
galaxies (Kravtsov 2013). The forced quenching simulations 
(red filled symbols) have 1.8 < 7 ' < 2 . 1 , while the standard 
simulations (black open symbols) have 2.3 < 7 ' < 2.4. Ob- 
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Figure 8. SUNRISE images of forced quenching simulations of halo4 at redshift 2 = 0 with resolution increasing from left to right. 


served galaxies exist with all such values, though the forced 
quenching simulations have more typical values, and while 
the forced quenching simulations fall within the observed 
distribution of y' and E s tar, the standard simulations fall 
well outside. 

It is also possible to measure the mass slope over a 
longer baseline using the velocities at r 1 / 2 and Raoo- Obser¬ 
vations using weak gravitational lensing and satellite kine¬ 
matics find that for massive ellipticals Vj / 2 ~ V 200 (Dut¬ 
ton et al. 2010). The standard simulations have declining 
circular velocity profiles with V1/2/V200 ~ 1.6 to 2.0. The 
forced quenching simulations have much flatter profiles with 
V 1 / 2 /V 200 ~ 0.8 to 1.0 in agreement with observations. 

The slope of the mass profile both locally at ~ rq / 2 and 
globally between rq/ 2 and R 200 further support the conclu¬ 
sion that the FQ simulations have realistic distributions of 
stars and total matter at small radii, while the standard 
simulations have too much stellar and total mass at small 
radii. 


3.4 Resolution effects on galaxy structure 

One of our initial conditions (halo4) has been run at four 
different resolution levels: levelO is the lowest and level3 is 
the highest (see Fig. |T] & Table [2]). In Figs. [4][T] the galaxy 
structural parameters for these simulations are shown with 
larger circles for higher resolution. There is clearly signifi¬ 
cant variation in the galaxy structural parameters at differ¬ 
ent resolutions. Here we focus our discussion on the forced 
quenching simulations at 2 = 0. 

The stellar mass increases with resolution by 0.39 dex 
from lowest to highest. However, most of this difference (0.33 
dex) is due to the lowest resolution simulation, which has 
only ~ 10 5 dark matter particles within the virial radius. 
There is just 0.06 dex difference between levell and level3. 
Circular velocity at the half-mass radius, Vi/ 2 , shows a simi¬ 
lar qualitative trend, with higher velocities in higher resolu¬ 
tion simulations, but with smaller relative changes. Between 
the highest three resolution simulations the variation in V 1/2 
is just 2%. Half-mass sizes vary by a factor of ~ 2 with no 
clear trend with resolution. Stellar surface densities show 
the largest variation with a factor of ~ 8 difference. This 
variation in surface density is clearly visible in the galaxy 
images shown in Fig. [8] If we remove the lowest resolution 
simulation the variations in sizes and surface densities are 
reduced to a factor of ~ 1.6 and ~ 2.0, respectively. 

I 11 summary, for simulations with more than a million 
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Figure 9. Circular velocity profiles for most massive galaxies in 
the standard simulations at 2 = 2.1 (solid lines): Total (black); 
dark matter (red); baryons (blue); stars (cyan); gas (green). The 
dashed lines show the dark and baryons for the no cooling simu¬ 
lation, which acts as a control to the impact of galaxy formation. 
The filled symbol corresponds to the circular velocity at the stel¬ 
lar half-mass radius. 


particles inside the virial radius the stellar masses and circu¬ 
lar velocities are well converged while the half-mass sizes and 
average stellar surface densities show significant variation. 


4 DARK HALO RESPONSE TO GALAXY 
FORMATION 

In the previous section we have established how well our 
various simulations reproduce observed global properties of 
elliptical galaxies at 2 = 0 and their progenitors at 2 ~ 2. 
We now turn our attention to the radial distribution of dark 
matter, and in particular how this responds to the galaxy 
formation process. 
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Figure 10. Circular velocity profiles for most massive galaxies 
in the standard simulations at z = 0 (solid lines): Total (black); 
dark matter (red); baryons (blue); stars (cyan); gas (green). The 
dashed lines show the dark and baryons for the no cooling simu¬ 
lation, which acts as a control to the impact of galaxy formation. 
The filled symbol corresponds to the circular velocity at the stel¬ 
lar half-mass radius. 


Figure 11. Circular velocity profiles for most massive galaxies 
in the forced quenching simulations at z = 0 (solid lines): To¬ 
tal (black); dark matter (red); baryons (blue); stars (cyan); gas 
(green). The dashed lines show the dark and baryons for the no 
cooling simulation, which acts as a control to the impact of galaxy 
formation. The filled symbol corresponds to the circular velocity 
at the stellar half-mass radius. 


4.1 Circular velocity profiles 

Fi gs. EpTI s how circular velocity profiles, where Wire = 
^ /GM{r)/r , for our simulated galaxies at z = 2.1 and z = 0. 
Velocity profiles are plotted from r m i n to the the virial ra¬ 
dius, R, 2 qo- ffere r m in is the maximum of the convergence 
radius as defined by Power et al. (2003) and twice the soft¬ 
ening length. The circular velocity at stellar half-mass radius 
is shown with a different symbol for each halo, following on 
from previous plots. 

When making these velocity profiles we have taken care 
to ensure that the systems are relaxed. During a merger 
event the assumption of spherical symmetry breaks down, 
the center is no longer well defined, and the derived mass 
profile will have an apparent core in the center. By compar¬ 
ing the mass profiles at different snapshots spaced nearby 
in time, the presence of mergers is easily detected. We have 
selected the snapshot closest to z = 0 or 2 = 2.1 in which 
both the galaxy formation and no cooling simulations show 
no obvious signs of being compromised by mergers. 

A comparison of the circular velocities at the virial ra¬ 
dius (the last point in the velocity profile) shows that the 
standard (solid lines, Fig. 1101) . forced quenching (solid lines, 
Fig. 111!) and no cooling (dashed lines, Figs. llOl fc llll) simula¬ 
tions all have the same dark and baryonic masses. Further¬ 
more, this implies that the haloes have been able to retain 
all of their cosmic share of baryons. In all simulations stars 
(cyan lines) dominate the baryon budget at small radii, while 
gas (green lines) becomes as important as the stars at large 
radii. In the FQ simulations the gas roughly follows the mass 



Figure 12. Change in dark matter circular velocity due to galaxy 
formation: with standard physics at 2 = 2.1 (blue, short-dashed 
lines) at z = 0 (black, solid lines); and with forced quenching at 
z = 0 (FQ, red, long-dashed lines). The filled symbols indicate 
the half-mass radii of the stars. All simulations contract within 
~ 20 kpc, but there is substantially less contraction in the FQ 
simulations. 
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Figure 13. Change in enclosed dark matter mass due to galaxy formation: with standard physics at z = 2.1 (left) at z = 0 (middle); and 
with forced quenching at z = 0 (right). Predictions from the adiabatic contraction (“AC”) models of Blumenthal et al. (1986, B86) and 
Gnedin et al. (2004, G04) are shown with green dotted and magenta dot-dashed lines, respectively. At redshift z = 0 the contraction in 
our simulations is substantially weaker than predicted by these models. Upper panels show change with respect to simulations without 
cooling (“No Cool”), while lower panels show differences with respect to the galaxy formation simulations (“Hydro”). 


profile of the dark matter, but in the standard simulations 
there is an additional component of concentrated cold gas. 

A comparison between the solid and dashed red lines 
in Figs. [9]- QT] shows that all the dark matter haloes have 
contracted in response to galaxy formation. The strength of 
the contraction is more clearly seen in Fig. [T^] which shows 
the ratio of dark matter circular velocities. The numerator 
is the dark matter circular velocity in the simulations with 
galaxy formation (i.e., standard/FQ) and the denominator is 
the dark matter circular velocity in the simulations without 
galaxy formation (i.e., no cooling). In the standard simula¬ 


tions the velocity ratio is roughly the same at z = 0 (black 
lines) as it is at z = 2.1 (blue lines), whereas in the FQ 
simulations (red lines) the contraction is not as strong. 

4.2 Adiabatic contraction formalism 

To understand the effect of the baryonic physics on the dark 
matter density profile, we employ the analytic adiabatic con¬ 
traction formalism outlined in Blumenthal et al. (1986), and 
introduced in another context in Barnes & White (1984). 
The main assumption is that the time scale for galaxy for- 
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mation is long compared to the orbital time scale of the dark 
matter particles. With the simplifying assumptions of spher¬ 
ical symmetry and circular dark matter particle orbits, the 
adiabatic invariant reduces to rM(r) = const , where M(r) 
is the total mass enclosed within radius r. 

Thus given a spherically enclosed mass profile from a 
simulation without gas cooling, Mi(r), (where the i refers 
to initial), we can derive the final dark matter profile, 
A/dm,f(r), once the final baryonic mass profile, Mb ar ,f(r), 
is specified. The initial total mass profile is split into dark 
matter Md m ,i(r), and baryons Mb ar ,i(r) either explicitly 
(as in our no cooling simulations) or implicitly assuming 
-T/bar, i (r ) = (1 — /bar) Aft (r) , where /bar is the cosmic baryon 
fraction (flb/flm ~ 0.16). 

An additional assumption is that the dark matter shells 
do not cross: M dm , i( r t) = Md m ,i(ri), where n is the “ini¬ 
tial” radius of a shell of dark matter, and rf is the radius of 
this shell after the effects of galaxy formation are included. 
Putting this together yields 

rt/n = Mi(n)/[M b ,f(r f ) + M d ,i(n)]. (3) 

Thus given Aft and Mb,f, one can solve Eq.[3]for the mapping 
between rf and n, and hence derive the final dark matter 
profile. The appeal of this formalism is that in the adiabatic 
limit the response of the halo depends only on the final state 
of the baryons, and is independent of how it was assembled. 

Particle orbits in CDM haloes are not circular. To ac¬ 
count for this Gnedin et al. (2004) introduced a modified adi¬ 
abatic invariant: rM(f ), where r and r are the current and 
orbit-averaged particle positions. The orbit average radius 
can be approximated as: f « f? v i r A(r/ R v n) w , with A ~ 0.85 
and w « 0.8. This modified adiabatic invariant results in 
slightly weaker contraction. 

Fig. [13] shows the change in the dark matter masses 
profiles in our simulations and that predicted by the adia¬ 
batic contraction models of Blumenthal et al. (1986, B 86 ) 
and Gnedin et al. (2004, G04). At z = 2.1 (left panel) the 
G04 model works for radii r > 5 kpc, and the B 86 model for 
r > 10 kpc. At smaller radii the models predict more con¬ 
traction than seen in our simulations. By 2 = 0 (middle and 
right panel) the over prediction has grown and extends to all 
radii (less then the virial radius). For example, at 5 kpc the 
G04 and B 86 models over predict the mass by a factor of ~ 2 
and ~ 3, respectively. These results are qualitatively simi¬ 
lar to previous studies on Milky Way mass haloes (Abadi 
et al. 2010; Pedrosa et al. 2010). 

Returning to the halo response formalism, an example 
of how the radii and masses are calculated in one of our 
simulations is shown in Fig. 1141 The solid lines show the 
total (black) and dark matter (red) mass profiles for the 
standard simulation for halo4.2. The dashed lines show the 
total (black) and dark matter (red) mass profiles for the 
corresponding simulations with no cooling. The horizontal 
dotted blue line shows an arbitrary mass (here 10 11 Mq), 
which intersects the final dark matter mass profile at rf and 
the initial dark matter mass profile at n. Mf and Mi are 
then the total enclosed mass at rf and n , respectively. This 
process is repeated for different arbitrary masses to obtain 
the relation between rf /n and Mi/Mf. 

Fig. 1151 shows the tracks of our simulations in the n/n 
vs Mi/Mf plane. The y-axis shows the “contraction” factor, 
while the x-axis can (typically) be mapped monotonically to 



Figure 14. Example showing how rf,r;,Mf,M; are calculated 
in the cumulative mass vs radius plot for simulation Halo4 using 
standard galaxy formation physics. The horizontal blue dotted 
line shows an arbitrary mass, here 10 11 Mq . The radius this in¬ 
tersects the dark matter profile from the standard simulation is 
termed rf (vertical blue solid line), while the radius this intersects 
the dark matter profile from the no cooling simulation simulation 
is termed n (vertical blue dashed line). The total masses con¬ 
tained within rf and rf are then given by Mf (horizontal blue 
solid line) and Mi (horizontal blue dashed line) respectively. 

radius. The adiabatic contraction formula (Eq. [3]) predicts 
that rf/n = Mi /Mf , which is indicated by the diagonal dot¬ 
ted line in the figure. No change in the dark matter profile 
corresponds to the horizontal dotted line at rf/n = 1, and 
expansion to rt/n > 1. 

All simulations show contraction at small radii (see also 
Figs. HUB, with the standard simulations at z = 0 show¬ 
ing a small amount of expansion at large radii, r > 30 
kpc, which corresponds to the radius where there is equal 
baryons and dark matter. The points show the halo response 
at rt = rf/ 2 . The standard simulations at 2 = 2.1 and 2 = 0 
have rf/n ~ 0.5, vs ~ 0.1 for adiabatic contraction, and 
the FQ simulations have rt/n ~ 0.8 vs ~ 0.4 for adiabatic 
contraction. 

Overall the halo response does not follow a single track, 
although interestingly at a given redshift each type of simu¬ 
lation results in a similar halo response. To quantify different 
halo responses we first consider the equation introduced by 
Dutton et al. (2007): 

rt/n = (Mi/Mf) 1 ' . (4) 

Adiabatic contraction is thus v = 1, contraction 0 < v < 1, 
no change to v = 0, and expansion v < 0. The cases v = 
0, 0.25, 0.5 and 1 are shown in Fig. [15] with dotted green 
lines. We note that the model of Gnedin et al. (2004) can 
be approximated with v ~ 0.8, and thus can be considered 
moderately strong contraction. 
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Figure 15. Dark halo response of our simulations: standard 2 = 2.1 (upper left, blue); standard 2 = 0 (upper right, black); forced 
quenching 2 = 0 (lower right, red); forced quenching 2 = 0 resolution tests with halo4 (lower left, red). The filled points show the halo 
response at the galaxy half-stellar mass radius (with different simulations represented by different symbols as indicated). The dotted green 
lines show various halo responses of the form n/n = (Mj/Mf) 1 '. Where v = I corresponds to adiabatic contraction, v = 0 corresponds 
to no change, and n = 0.5 and n = 0.25 correspond to weak contraction. The dot-dashed magenta lines show the halo response model of 
Abadi et al. (2010, A10). 


While the standard simulations at 2 = 0 are approx¬ 
imated by Eq. [4] with v = 0.25, the standard simulations 
at 2 = 2.1 and the FQ simulations at z = 0 are not well 
described by this formula. The magenta lines show the for¬ 
mula from Abadi et al. (2010), rt/n = 1 — 0.3(Mi/Mf — l) 2 , 
which approximately describes our forced quenching simu¬ 
lations at 2 = 0 (lower right panel). This similarity may be 
a coincidence as their simulations are missing important as¬ 
pects of galaxy formation processes such as star formation 
and feedback. 

The shaded regions in Fig. [15] bracket the halo response 
shown in the 4 initial conditions. They are given by the fol¬ 


lowing relations, where the uncertainties bracket the range 
in halo response: standard simulations at 2 = 2.1 (upper 
left) 

rt/n = 0.5(i°;}° 2 ) + 0.5(t°- 2 )(Mi/M f + O.Ol/; 2 ) 2 ; (5) 

standard simulations at 2 = 0 (upper right) 
rt/n = 1.0(1°;°®) - 0.52(Mi/M f - l) 2 ; (6) 

FQ simulations (lower right) 

rt/n = 0.75(1°;°7) + 0.25(;°03)(M/M f ) 2 ; (7) 

A numerical convergence test is shown in the lower left 
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Figure 16. Logarithmic slope of the dark matter density profile, 
a, measured between 1 and 2% of the virial radius, f?200, vs the 
ratio between the stellar and virial mass, M st ar/AL200- Blue sym¬ 
bols show standard simulations of halo4 at z = 2.1, black symbols 
show standard simulations at z = 0, and red symbols show forced 
quenching simulations at z = 0. The green dotted lines show 
a ~ —1.2 for the control simulations (without cooling). The grey 
shaded region shows the relation from Toilet et al. (2015, in prep) 
which was calibrated against haloes of mass 10 10 — 10 12 Mq. 

panel of Fig. 1151 This shows the halo response for halo4 run 
at four different mass (and force) resolution levels (see Ta¬ 
ble [2J. All simulations have similar halo response (especially 
the highest three resolution simulations), suggesting conver¬ 
gence. The lowest resolution simulation (halo4.0) does not 
have as much contraction, but follows a similar path in the 
n/r\ vs Mi/Mt plane. 

4.3 Central dark matter slopes 

Another way to express the strength of the contraction of 
the dark matter halo is the central density slope, a. We 
measure a between 0.01 and 0 . 02 f ?200 which corresponds to 
~ 5 — 10 kpc for our simulations at 2 = 0. Following Di 
Cintio et al. (2014) we show in Fig. [16] the relation between 
a and galaxy formation efficiency, M B t al /M 2 oo- Our FQ sim¬ 
ulations (red points) have —1.7 < a < —1.4, while the stan¬ 
dard simulations (black points) have —2.0 < a < —1.6. For 
reference, the dotted green lines show the no cooling sim¬ 
ulations which have —1.3 < a < —1.1, the same values as 
found in dark matter only simulations. 

The solid line and shaded region shows the relation ob¬ 
tained from Toilet et al. (2015, in prep) who used a set of 
~ 80 zoom-in simulations of halo mass ~ 10 10 to ~ 10 12 from 
the NIHAO project (Wang et al. 2015). Our standard simu¬ 
lations are close to this relation (both at 2 = 0 and 2 = 2.1), 
but the forced quenching simulations are significantly offset. 
At a stellar to halo mass ratio of Af sta r/A /200 ~ 0.01 spiral 


galaxy haloes have large amounts of expansion, whereas our 
elliptical simulations contract. Thus the inner dark matter 
density slope is not purely determined by the efficiency of 
star formation. This difference is a reflection of the different 
formation channels of spiral vs elliptical galaxies. In spiral 
galaxies the vast majority of stars form in situ, with the effi¬ 
ciency being regulated by SN feedback which also has strong 
expansive effects on the halo. In elliptical galaxies, a signif¬ 
icant fraction of the stars were formed ex situ, and were 
assembled with dry mergers which has a milder expansive 
effect on the halo. 

Finally, we note that the form of the halo response in 
Eqs.[5]0 which has a constant ratio between initial and final 
radii at r = 0, implies that at very small radii the central 
dark matter slope tends to the dissipationless case. If haloes 
are described by the Navarro, Frenk & White (1997) for¬ 
mula, then the inner slope would be —1, however, if the 
Einasto (1965) profile holds then this implies the centers of 
dark matter haloes have constant density. Practically speak¬ 
ing, however, the presence of a supermassive black hole may 
cause a dark matter spike (Gondolo & Silk 1999), rendering 
these extrapolations a moot point. 

4.4 Slow expansion 

When we disable gas cooling at 2 = 2.1, the gas at small radii 
heats up and expands to larger radii, see the green lines in 
Fig. 1171 One might worry that this would cause the halo 
to expand artificially. However, the circular velocity plots in 
Fig- E] show that the gas is a sub-dominant component with 
Vgas ~ 0.1 Vcirc, and thus contributes just ~ 1% to the en¬ 
closed mass at all radii (recall that M oc V 2 ). Thus in the ex¬ 
treme case that all of the gas is removed instantaneously, the 
expansion of the dark and stellar mass distributions would 
be negligible. 

Fig. 1171 shows the evolution in the mass enclosed within 
5 proper kpc for our simulations. The total mass is shown in 
black, the dark matter in red, the stars in blue and the gas 
in green. The three line types correspond to the standard 
(solid), FQ (long-dashed), no cooling (dotted) simulations. 
Time steps where the halo is clearly undergoing a merger 
have been removed since it is impossible to determine the 
central 5 kpc accurately during those time periods. 

I 11 the no cooling simulation the dark matter mass is 
roughly constant at ~ 2 —3 x 10 10 Mq since a = 0.15 (2 ~ 5). 
In the standard simulation the dark halo starts to contract 
at 2 ~ 5, reaches a maximum mass a factor of ~ 2.5 times 
that of the no cooling simulation at 2 ~ 2, and remains 
roughly constant to 2 = 0. By contrast the stellar mass 
within 5 kpc continues to increase all the way to 2 = 0. In 
the FQ simulation the stellar mass is roughly constant since 
2 = 2.1, with a slight reduction due to stellar mass loss, 
while the dark matter mass slowly decreases since 2 ~ 2 by 
~ 0.2 dex. This gradual expansion of the dark matter halo 
since 2 = 2.1 in the forced quenching simulations is further 
evidence that the expansion is not caused by a single event 
at 2 = 2.1 or later, such as a major merger. 

There are two processes that are likely contributing to 
the smooth expansion: multiple minor mergers (e.g., El-Zant 
et al. 2001) and adiabatic expansion due to stellar mass 
loss. In our simulations up to 40% of the stellar mass, i.e., 
~ 4x 10 1 o Mq, is returned to the ISM. I 11 the standard simu- 
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Scale Factor a = 1/(1+z) Scale Factor a = 1/(1+z) 


Figure 17. Mass enclosed within 5 proper kpc vs scale factor for the most massive progenitors in our 4 simulations. The total mass 
is given by black lines and points, dark matter by red and stars by blue. In the no cooling simulations (crosses, dotted line) the dark 
matter is in place by a ~ 0.15 [z ~ 5). In the standard simulations (filled circles, solid lines) the halo contracts, increasing the mass by 
a factor ~ 2.5. The contraction occurs by z = 2.1 and remains constant to z ~ 0, even though the stellar mass increases by an order of 
magnitude during thus interval. In the forced quenching simulations (open circles, dashed lines) the dark matter within 5 kpc gradually 
decreases with time due to a combination of stellar mass loss mass and dry minor merging. 


lations this gas gets turned into new stars, but in our forced 
quenching simulations this gas gets heated and expands to 
larger radii resulting in adiabatic expansion. This stellar re¬ 
cycling explains why the stellar masses of our forced quench¬ 
ing simulations declines gradually with time, except for the 
occasionally merger. The signature of the adiabatic expan¬ 
sion processes is that within a fixed radius the relative re¬ 
duction in dark matter should equal the relative reduction 
in stellar mass. This appears to be the case for halol and 
halo2, but halo3 has a merger event at a ~ 0.45 which in¬ 
creases the stellar mass, but decreases the dark mass, and 
halo4 has stronger reduction in dark mass than stellar mass. 


We thus conclude that both minor merging and stellar 
mass loss are responsible for the expansion of the dark mat¬ 
ter halo since z = 2.1 in our forced quenching simulations. 


5 DISCUSSION 

For the massive galaxies we simulate in this paper there are 
three main processes that determine the halo response: dissi¬ 
pative (gas) accretion; dissipationless (dry) merging and gas 
outflows driven by AGN. In lower mass galaxies SN/stellar 
winds are the dominant drivers of gas outflows. The first 
two processes are included in our simulations, while AGN 
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are not. Dissipative accretion will only cause the halo to con¬ 
tract, while dry merging and AGN feedback can (but will not 
necessarily) cause halo expansion. These processes are not 
independent, as for example, stronger dissipation leads to 
more contraction in the central galaxy, but also denser satel¬ 
lites which can survive to smaller radii, and cause more ex¬ 
pansion to the dark matter. Thus, it is not obvious whether 
more dissipation will always translate into more contraction. 

5.1 Does dry merging lead to halo expansion? 

In the halo quenching scenario the mass assembly is dom¬ 
inated since z ~ 2 by dry mergers. It has been proposed 
(El-Zant et al. 2001; Lackner & Ostriker 2010) that this 
will create cores in the dark matter halo. While we do find 
that dry merging causes haloes to expand relative to the 
z ~ 2 case, the net effect at the halo mass scale we study 
( M 200 ~ 1O 13 M 0 ) is still contraction: the dark mass within 
5 kpc has increased by factor of ~ 1.4. Our simulations show 
that the contractive effects of dissipation at early times out¬ 
weigh the expansive effects of dry merging at late times. In 
contrast to previous studies, our simulations are fully cos¬ 
mological with realistic progenitor masses and sizes which 
form realistic elliptical galaxies at z = 0. This is important, 
as the stellar density, orbits, and number of satellite galaxies 
determines the magnitude of the halo expansion effect. 


5.2 What is the halo response due to AGN 
feedback? 

Simulations that incorporate AGN driven feedback can re¬ 
sult in halo expansion (e.g., Peirani et al. 2008; Duffy 
et al. 2010; Martizzi et al. 2012), although others seems to 
result in no significant change (e.g., Schaller et al. 2015a). 
Since our forced quenching simulations produce contracted 
haloes, this difference in halo response suggests a new way to 
distinguish between suppressive (halo) quenching and ejec- 
tive (AGN) quenching. 

It still needs to be determined what simulations with 
AGN feedback actually predict for halo structure. It is im¬ 
portant that the simulations reproduce the global proper¬ 
ties of galaxies such as we discuss in § [3] since the strength 
of the feedback will effect both the halo response and the 
masses and structure of the galaxies. Theoretical models of 
AGN feedback are often calibrated to match the present 
day galaxy stellar mass function (e.g., Schaye et al. 2015). 
A complication to this calibration is that the true stellar 
masses of galaxies are not known accurately. In the case of 
the EAGLES simulations (Schaye et al. 2015) the observed 
stellar mass function used in the calibration was derived 
assuming a Milky Way IMF. We know from dynamics and 
lensing studies that a Milky Way IMF requires close to adia- 
batically contracted dark matter haloes in elliptical galaxies 
(Auger et al. 2010a; Schulz et al. 2010; Dutton et al. 2011b; 
2013b). However, the halos in the EAGLES simulations show 
almost no change (Schaller et al. 2015a), which points to¬ 
wards an inconsistency between these simulations and ob¬ 
servations. 

Furthermore, recent evidence points towards non-Milky 
Way IMFs in massive elliptical galaxies from both dy¬ 
namical and stellar population modeling approaches (e.g., 


Auger et al. 2010a; van Dokkum & Conroy 2010; Conroy 
& van Dokkum 2012; Dutton et al. 2013a,b; Cappellari 
et al. 2013; Barnabe et al. 2013; Ferreras et al. 2013; Sonnen- 
feld et al. 2015). The favored stellar masses are a factor of 
~ 2 higher than obtained assuming a Milky Way IMF. Thus 
it would seem necessary for simulations such as EAGLES to 
be re-calibrated to galaxy mass functions based on heavier 
IMFs in the most massive galaxies. 

The structural properties of galaxies can provide addi¬ 
tional constraints for AGN feedback models. In particular 
the total mass density slopes at ~ 1% of the virial radius, 
and the size vs velocity relation are independent of uncer¬ 
tainties in the IMF. At a halo mass of ~ 3 x 1O 13 M 0 the 
total mass density slopes at 1% of the virial radius in the 
EAGLES simulations are typically shallower than observed 
(Schaller et al. 2015b), providing further evidence for insuf¬ 
ficient dissipation. 

Simultaneously reproducing several galaxy structural 
properties is a non-trivial task. Dubois et al. (2013) present 
six zoom-in simulations of halo mass 0.4 — 8 x 1O 13 M 0 that 
include AGN feedback. The simulations broadly reproduce 
the size vs stellar mass, velocity dispersion vs stellar mass 
and total mass density slopes at redshift a ~ 0. In detail 
there are some discrepancies that again point towards in¬ 
sufficient dissipation. The total mass density slopes are too 
shallow, and relative to observed stellar masses derived us¬ 
ing a Salpeter IMF the simulated sizes are too large while 
the velocity dispersions are slightly too small. 

5.3 Future directions for halo quenching models 

The halo quenching model we implement in this paper could 
be improved upon in several ways. We turn the gas adiabatic 
everywhere in the simulation at z ~ 2, which corresponds 
to when the most massive progenitor reaches a halo mass 
of ~ 1O 12 M 0 . A more realistic scenario would be to only 
shut off cooling in haloes above ~ 1O 12 M 0 , or some other 
galaxy property such as stellar density or bulge mass, al¬ 
lowing cooling and star forming in lower mass progenitor 
galaxies to occur below z = 2.1. This additional cooling 
would presumably result in denser inner galaxies, and more 
contraction. Another improvement would be to implement 
photo-ionization from stellar sources (Kannan et al. 2014, 
2015) or heating from AGB stars (Conroy et al. 2015) di¬ 
rectly into the simulations. 

Simulations could be run at different halo masses, e.g., 
from ~ 10 12 — 1O 14 M 0 to determine if there is any mass 
dependence to the halo contraction in halo quenching sce¬ 
nario. In the context of ACDM, more massive haloes have 
a larger fraction of the stellar mass in the central galaxy 
assembled through mergers (Behroozi et al. 2013), and thus 
the expansive effects of dissipationless assembly would be 
expected to be more important in higher mass haloes. In¬ 
deed, at the galaxy cluster scale M 200 ~ 1O 15 M 0 Laporte & 
White (2015) find that dissipationless assembly since z = 2 
reverses the contraction in the progenitor galaxies resulting 
in net halo expansion within ~ 1% of the virial radius by 
redshift z = 0. We note that at a halo mass of 1O 13 M 0 
the accreted fraction of stars is already quite high at ~ 60% 
(Behroozi et al. 2013), thus any mass dependence is likely to 
be weak. Furthermore, since the halo and galaxy mass func¬ 
tions are steep, very few quiescent galaxies live in the highest 
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Figure 18. Dark matter fraction vs circular velocity, both mea¬ 
sured at 5 kpc. Symbols show our simulations while lines show 
observed values for different halo response models from Dutton 
et al. (2013b): no contraction (v = 0), weak contraction (y = 0.5) 
adiabatic contraction (AC, v = 1). The forced quenching simula¬ 
tions (red) are most similar to the weak contraction model. The 
standard simulations (black) have very low dark matter fractions, 
which observationally requires uncontracted halos, but the simu¬ 
lations have contracted indicating an inconsistency. 


mass dark matter haloes, and thus majority of massive el¬ 
liptical galaxies will have contracted dark matter haloes in 
the halo quenching scenario. 


5.4 Observations of dark matter fractions 

An observational test of halo response models is the dark 
matter fraction on scales of galaxy half-light radii, i.e., ~ 5 
kpc. Total masses on these scales can be measured reliably 
using stellar kinematics and/or strong gravitational lensing 
(e.g., see Courteau et al. 2014 for a review). However, de¬ 
composing the total mass into baryons and dark matter is 
more challenging. Stellar population synthesis models can be 
used to convert luminosity profiles into stellar mass profiles 
up to the stellar IMF and systematics uncertain phases in 
stellar evolution. It is well established observationally that 
if the IMF is universal, then massive elliptical galaxies have 
dark matter fractions of /dm {Re) ~ 0.5 and in the context 
of ACDM, require close to adiabatic halo contraction (Auger 
et al. 2010a; Schulz et al. 2010; Dutton et al. 2011b, 2013b). 

In Fig. [TS] we show the dark matter fractions within 
spheres of radius 5 kpc vs the circular velocity at 5 kpc 
in our simulations (symbols) and compared to observation- 
ally constrained models (lines) from Dutton et al. (2013b). 
This plot is similar to Fig. 21 of Dubois et al. (2013), but 
we use a fixed physical radius, rather than a relative radius 
such as the half-light radius to eliminate any differences of 
scale between the simulations and observations. The forced 


quenching simulations have dark matter fractions of ~ 30%, 
while the standard simulations have ~ 10%. The lines show 
models with different assumptions about the halo response 
parameterized by Eq. [4] from no change (v = 0) to adia¬ 
batic contraction (u = 1). Models with stronger contraction 
have correspondingly lower stellar mass to light ratios. At 
a circular velocity of 300 km s^ 1 the stellar mass offsets are 
0.23, 0.16, and 0.03 dex, where 0.23 corresponds to a Salpeter 
IMF, and 0.0 to a Chabrier IMF. Note that all models re¬ 
produce (by construction) the scaling relations between size, 
stellar mass, and velocity dispersion. 

An accurate measurement of the dark matter within 5 
kpc would thus be able to falsify our simulations for elliptical 
galaxy formation. We note that variety of studies seem to 
favor models closer to v = 0 and Salpeter IMFs than v = 1 
and Chabrier IMFs (Auger et al. 2010a; Dutton et al. 2013b; 
Dutton & Treu 2014; Sonnenfeld et al. 2015), although there 
are some exceptions (Smith et al. 2015). Clearly more work 
needs to be done to make the distinction more definitive. 


5.5 Priors on dark matter density slopes 

Our results have implications for mass modeling of ellipti¬ 
cal galaxies. Observations commonly parameterize the dark 
matter halo as either a power-law with a free slope, or a 
double power-law with a free inner slope: 

p(r)/po = (r/r s )~ a (l + r/r s )~ 3+a . (8) 

For example, Cappellari et al. (2013) adopt Eq. [H] with a 
uniform prior on inner slope with limits [0,1.6]. The upper 
limit is consistent with our forced quenching simulations, 
but excludes simulations with stronger contraction. Due to 
the degeneracy between the halo response and dark mat¬ 
ter fraction with the stellar mass-to-light ratio (e.g., Dutton 
et al. 2013b) the choice of prior is important if one wishes to 
accurately constrain stellar mass-to-light ratios and make 
inferences on the IMF. Our simulations suggest an upper 
limit of a = 1.6 is too restrictive. A more conservative up¬ 
per limit of a = 2 allows for the full range of inner slopes 
found in our simulations. 


6 SUMMARY 

We use cosmological hydrodynamical simulations with the 
SPH code GASOLINE to investigate the response of dark mat¬ 
ter haloes to the formation of massive elliptical galaxies. We 
consider 4 initial conditions that grow into haloes of mass 
M 200 ~ 10 13 Mq by the present day. At our standard resolu¬ 
tion each simulation has more than 1.5 million dark matter 
particles within the virial radius at a = 0 enabling us to 
resolve the dynamics at 1% of the virial radius. Our highest 
resolution simulation has 10.6 million dark matter particles 
within the virial radius at 2 = 0, making it one of the highest 
resolution simulations of elliptical galaxy formation to date. 
Our standard simulations include metallicity dependent gas 
cooling, star formation, and stellar feedback. 

We summarize our results as follows: 

• At redshift 2 = 2.1 our standard simulations have 
stellar to halo mass ratios consistent with halo abundance 
matching, assuming a Salpeter IMF to derive the stellar 
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masses (Fig. [3}. Galaxy half-mass sizes and circular veloci¬ 
ties are also consistent with observations (Fig. [4j). 

• At 2 = 2.1 the dark matter haloes have contracted in 
response to galaxy formation (Figs. [9] 1121 1151) . but not as 
strong as predicted by the adiabatic contraction formalism 
(Blumenthal et al. 1986). 

• By z = 0 the standard simulations have overcooled, 
with a factor of > 4 times too many stars, resulting in cir¬ 
cular velocities, stellar densities, and mass density slopes 
that are too high (Figs. [3]-[7]) 

• We investigate the halo quenching scenario by shutting 
down cooling and star formation at 0 = 2.1 (when the most 
massive progenitors have M200 ~ 10 12 Mq) and evolving the 
simulation hierarchically to z = 0. We refer to these simula¬ 
tions as forced quenching (FQ). The resulting galaxies have 
many properties consistent with observed elliptical galaxies: 
M star ~2x 1O 11 M 0 , M s tar/M 2 oo ~ 0.01 (Fig. [3]), flat circu¬ 
lar velocity profiles with mass density slope 7 ' ~ 2 (Figs. [S] 
& 1111) . half-mass sizes of rq / 2 ~ 4 — 10 kpc (Fig. [4J, and 
circular velocities Kirc(ri/ 2 ) ~ 300 — 400kms -1 (Fig. H}. 

• In all the FQ simulations the dark matter haloes 
contract, although less than in the standard simulations 
(Fig. 1121) . The contraction is much weaker than predicted by 
the adiabatic contraction models of Blumenthal et al. (1986) 
and Gnedin et al. (2004), but can be described with a simple 
formula (Eq. [7]). 

• The dark matter density slopes (at redshift z = 0) mea¬ 
sured between 1 - 2 % of the virial radius vary from — 1 . 2 ± 0.1 
in the control simulations to —1.53 ± 0.17 in the FQ sim¬ 
ulations and to —1.76 ±0.17 in the standard simulations 
(Fig. [[ 6 ]). 

• Dry merging alone is unable to reverse the contractive 
effects of early dissipation, which must occur to form an 
elliptical galaxy. 

• Simulations in the literature find that AGN feedback 
can cause halo expansion (e.g., Duffy et al. 2010; Martizzi 
et al. 2012 ), and thus there may be qualitatively differ¬ 
ent halo responses, and dark matter fractions within ~ 5 
kpc, in the suppressive (halo) quenching and ejective (AGN) 
quenching scenarios. We note that different halo responses 
will also require different stellar IMFs in order to be con¬ 
sistent with observed constraints on dark matter fractions 
(Fig. [T51) . 

While our treatment of the halo quenching process is 
clearly an oversimplification, we hope that the excellent 
agreement between the structural properties of our forced 
quenching simulations with observations motivates further 
studies of the halo quenching scenario. 
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APPENDIX A: GALAXY PARAMETERS 

Table Eli lists the parameters of the most massive galaxy in 
each simulation at redshift z ~ 2 and z ^ 0. Note that the 
actual redshifts vary slightly to avoid major mergers that 
occur at the nominal output redshifts. 
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Table Al. Galaxy parameters. Column (1) name of the initial conditions - see Table [2] for simulation parameters. Column (2) type 
of simulation: standard (S) or forced quenching (FQ). Column (3), z, redshift of output. Column (4), M 200 , virial mass. Column 
(5), Mgtar, stellar mass inside 20% of the virial radius. Column ( 6 ), R\/ 2 -> 3D half-mass radius. Column (7), V 1 / 2 , circular velocity 
(V(r) = yjGM(r)/r) at the 3D half-mass radius. Column ( 8 ), V 200 ? circular velocity at the virial radius. Column (9), V^ pc , circular 
velocity at 5 kpc. Column (10), /dm? dark matter fraction at 5 kpc. Column (11), a, slope of the dark matter density profile measured 
between 1 and 2% of the virial radius. Column (12), 7 ', slope of the total mass profile measured between 1 and 2% of the virial radius. 
Column (13), £ s tar, average surface density of the stars inside the 2D half-stellar mass radius. 


Name 

type 

z 

log 10 M 200 

log 10 Afstar 

Rl/2 

Vi /2 

V 200 

Vikpc 

/dm 

a 

l' 

log 10 Sstar 




[M 0 ] 

[M 0 ] 

[kpc] 

[kms -1 ] 

[kms -1 ] 

[kms -1 ] 




[M 0 pc -2 ] 

a) 

( 2 ) 

(3) 

(4) 

(5) 

( 6 ) 

(7) 

( 8 ) 

(9) 

( 10 ) 

(ii) 

( 12 ) 

(13) 

halol 

S 

0.12 

13.14 

11.80 

4.28 

599.3 

355.5 

589.5 

0.146 

-1.88 

2.37 

3.99 

halo 2 

S 

0.10 

13.14 

11.82 

4.68 

580.7 

353.7 

577.4 

0.130 

-1.93 

2.42 

3.93 

halo3 

s 

0.08 

13.14 

11.94 

5.74 

602.7 

353.5 

613.6 

0.110 

-1.58 

2.28 

3.88 

halo4.0 

s 

0.00 

13.37 

12.09 

4.21 

825.1 

414.2 

798.6 

0.094 

-1.55 

2.42 

4.29 

halo4.1 

s 

0.12 

13.33 

12.09 

5.08 

756.7 

409.6 

758.6 

0.108 

-1.61 

2.36 

4.13 

halo4.2 

s 

0.00 

13.34 

12.07 

5.13 

744.4 

403.9 

747.6 

0.135 

-1.67 

2.36 

4.10 

halol 

FQ 

0.00 

13.14 

11.32 

6.45 

343.3 

348.4 

326.0 

0.349 

-1.71 

1.80 

3.17 

halo 2 

FQ 

0.14 

13.13 

11.00 

3.59 

288.2 

354.2 

295.3 

0.344 

-1.63 

2.10 

3.36 

halo3 

FQ 

0.00 

13.19 

11.41 

7.95 

340.5 

362.3 

321.1 

0.288 

-1.36 

1.79 

3.09 

halo4.0 

FQ 

0.00 

13.36 

11.09 

13.79 

331.5 

410.9 

253.4 

0.564 

-1.51 

1.49 

2.26 

halo4.1 

FQ 

0.00 

13.35 

11.42 

7.01 

369.5 

408.0 

381.6 

0.310 

-1.43 

2.08 

3.18 

halo4.2 

FQ 

0.00 

13.33 

11.44 

7.78 

372.1 

401.2 

380.6 

0.295 

-1.52 

2.03 

3.11 

halo4.3 

FQ 

0.00 

13.43 

11.48 

11.08 

376.9 

432.3 

348.3 

0.363 

-1.57 

1.82 

2.84 

halol 

S 

2.22 

12.15 

11.01 

2.06 

384.7 

237.0 

398.0 

0.390 

- 1.20 

1.21 

3.83 

halo 2 

S 

2.08 

12.19 

10.92 

1.43 

404.6 

240.6 

373.8 

0.360 

-1.61 

1.58 

4.06 

halo3 

s 

2.08 

12.32 

10.94 

1.51 

385.1 

265.3 

363.8 

0.356 

-1.57 

1.63 

4.03 

halo4.0 

s 

2.07 

12.30 

10.79 

2.69 

320.4 

259.9 

343.3 

0.404 

-1.98 

0.99 

3.38 

halo4.1 

s 

2.07 

12.25 

11.14 

0.96 

582.2 

251.2 

435.0 

0.354 

- 1.86 

2.31 

4.62 

halo4.2 

s 

2.22 

12.22 

11.13 

1.30 

517.2 

251.6 

437.1 

0.383 

-1.76 

1.93 

4.35 

halo4.3 

s 

2.07 

12.42 

11.04 

1.94 

410.0 

286.8 

401.0 

0.402 

-1.59 

1.79 

3.92 
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